Dipole-active optical phonons in YTiC>3: ellipsometry study and lattice-dynamics 

calculations 



00 

o 
o 

> 
O 

(N 



a 

C 

o 
o 



> 
o 



oo 
o 



N. N. Kovaleva, 1 ' 2 A. V. Boris, 1 ' 2 L. Capogna, 3 - 4 J. L. Gavartin, 5 
P. Popovich, 1 P. Yordanov, 1 A. Maljuk, 6 A. M. Stoncham, 7 and B. Keimer 1 

1 Max- Planck- Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
Department of Physics, Loughborough University, Loughborough, LE11 3TU, UK 
3 CNR-INFM, CRS-SOFT and OGG Grenoble, 6 Rue J. Horowitz, BP 156 F-38042 Grenoble CEDEX 9, France 
* Institut Laue Langevin, 6 rue J. Horowitz BP 156 38042 Cedex 9 Grenoble, France 
5 Accelrys, Cambridge, CB4 OWN, UK 
6 Hahn-Meitner- Institute, Glienicker Str. 100, D-14109 Berlin, Germany 
7 London Center for Nanotechnology and Department of Physics and Astronomy, 
University College London, Gower Street, London WC1E 6BT, UK 
(Dated: November 26, 2008) 

The anisotropic complex dielectric response was accurately extracted from spectroscopic ellip- 
sometry measurements at phonon frequencies for the three principal crystallographic directions of 
an orthorhombic (Pbnm) YT1O3 single crystal. We identify all twenty five infrared-active phonon 
modes allowed by symmetry, 7Bi„, 9B2 U , and 9B3„, polarized along the o, b-, and a-axis, respec- 
tively. From a classical dispersion analysis of the complex dielectric functions e(w) and their inverses 
-l/e(w) we define the resonant frequencies, widths, and oscillator strengths of the transverse (TO) 
and longitudinal (LO) phonon modes. We calculate eigenfrequencies and eigenvectors of Bi„, B2 U , 
and B3 U normal modes and suggest assignments of the TO phonon modes observed in our ellip- 
sometry spectra by comparing their frequencies and oscillator strengths with those resulting from 
the present lattice-dynamics study. Based on these assignments, we estimate dynamical effective 
charges of the atoms in the YTiOa lattice. We find that, in general, the dynamical effective charges 
in YTiOs lattice are typical for a family of perovskite oxides. By contrast to a ferroelectric BaTi03, 
the dynamical effective charge of oxygen related to a displacement along the c-axis does not show 
the anomalously large value. At the same time, the dynamical effective charges of Y and ab-plane 
oxygen exhibit anisotropy, indicating strong hybridization along the a-axis. 

PACS numbers: 78.30.-j, 78.20.Bh, 61.82.Ms 



I. INTRODUCTION 



The richly diverse behavior of rare-earth titanates 
i?Ti03 with a single 3d electron in the nearly degen- 
erate Ti tig orbitals is a consequence of the interplay 
between electron correlations, spin-orbit coupling, and 
electron-lattice coupling including the Jahn Teller (JT) 
effect [H, H, i, 1, i, i 0, 1, H, M El Gl IH ■ Fine 
adjustment of the perovskite lattice in the titanates, in 
accordance with the principle of close packing, occurs 
mainly due to the cooperative rotations of the TiOe oc- 
tahcdra, which arc strongly dependent on the ionic radius 
of the rare-earth atom. The orbital structure that arises 
due to these distortions is determined by vibronic in- 
teractions between the lattice and the orbital subsystem. 
These octahedral distortions are also coupled to the mag- 
netic ground state in i?TiOa: upon decreasing the ionic 
radius of R from La to Y, the transition temperature of 
the G-type antiferromagnetic state observed in LaTiC>3 
(Tjv = 150 K) shifts to progressively lower temperatures, 
and the ground state ultimately becomes ferromagnetic 
with Tn = 30 K in YTi0 3 [U L^LOIJS ■ Point-charge 



P, [gl, [l5[ and density functional |l2l" 13 1 calculations 
produce many aspects of the physical properties of the 
titanates, including the magnetic ground states of the 
end-point compounds LaTiOa and YTiOs, based on the 



experimentally determined lattice structure. These cal- 
culations imply that the t^ g levels are split by ~ 100-200 
mcV, so that the Ti orbitals are almost frozen in the ex- 
perimentally relevant temperature range. However, some 
aspects of the magnetic [J, @] and orbital [HI excita- 
tion spectra are difficult to reproduce in the framework 
of these rigid-orbital models. This has stimulated theo- 
ries based on substantial fluctuations in the orbital sector 
[H, ||. While such fluctuations are not apparent in the 
specific heat (l9j . recent thermal conductivity measure- 
ments of i?Ti03 have been interpreted as evidence of or- 
bital fluctuations [2(| • Other recent studies of YTiOs sin- 
gle crystals have revealed anomalous temperature depen- 
dencies of the optical spectral weight [2l| and thermal ex- 
pansion [22| with an onset around 100 K, well above Tc, 
which may be manifestations of a temperature-dependent 
rearrangement of the tig orbitals. 

Because of the intimate relationship between lattice 
structure, orbital occupation, and magnetism, the lattice 
dynamics plays an important role in the unusual proper- 
ties of the rare-earth titanates. In particular, one would 
expect anharmonic behavior and damping of the most 
strongly coupled modes in the spin-lattice-orbital sector 
due to orbital fluctuations, as well as anomalous energy 
shifts of phonons around the magnetic transition temper- 
ature [23, Therefore experimental and theoretical 
studies of the lattice dynamics have the potential to be- 
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come quantitative diagnostics of the mechanisms under- 
lying the unusual phenomena in the rare-earth titanates. 
However, as little information about the phonon spectra 
and the assignment of the phonon modes of i?Ti03 is thus 
far available, this potential has not yet been exploited. 

The approach described here combines theoretical 
lattice-dynamics calculations with accurate spectroscopic 
ellipsometry measurements on well-characterized, nearly 
stoichiometric YTiC>3 single crystals of high quality. 
From the ellipsometry experiments we obtain the fre- 
quencies and oscillator strengths of all infrared (In- 
active TO phonon modes allowed by the Pbnm sym- 
metry. The lattice-dynamics calculations have been 
performed using the GULP code [2j|. We calculate 
the eigcnfrequcncics of the Brillouin-zone center normal 
modes of Bi u , B 2 „, and B3„ symmetry and the corre- 
sponding eigenvector components. Based on these calcu- 
lations, we suggest assignments of the TO phonon modes 
observed in our ellipsometry spectra by comparing their 
frequencies and oscillator strengths with the calculated 
values. In addition, we estimate the dynamical effec- 
tive charges of the atoms in the YTi03 lattice using the 
calculated eigenvector components and the experimental 
phonon parameters of the TO frequencies and oscillator 
strengths. The assignment is further verified by ab ini- 
tio lattice-dynamics calculations [26| within the density 
functional approach using the VASP software package 
[27l . [28^ . The present study provides a coherent picture 
of the lattice vibrations in YTi03, consistent with both 
experimental data and with full electronic structure cal- 
culations. 



II. CRYSTAL STRUCTURE AND T-POINT 
PHONONS IN YTi0 3 

YTi03 exhibits an orthorhombic structure of a 
GdFe03 type described by the Pbnm space group (Pnma 
D^h, No. 62 in the International Tables for X-ray Crys- 
tallography p9| . with nonconventional coordinate axes 
orientation, Pnma(a,b,c) Pbnm(b, c, a)), with 4 

f.u./unit cell (see Fig. 1, on the right side). This or- 
thorhombic structure can be regarded as a distorted mod- 
ification of the simple cubic ABO3 perovskite structure 
with space group Pm3m (see Fig. 1, on the left side). The 
orthorhombic distortions can be attributed to the com- 
mon atomic size mismatch: the sum of the ionic radii in 
the T1-O2 layer, m + ro, docs not match that of the 
Y-0 layer, (ry + ro)/V2, in the right way for a sta- 
ble cubic structure. These orthorhombic distortions lead 
first to the oxygen ion displacements, resulting from the 
two consequent and coordinated rotations of the TiOe 
octahedra, described by a~a~b + in the Glazer's nota- 
tions [3(|, and second to a modification of the Y-0 and 
Y-Ti coordinations. In addition, there are JT-type dis- 
tortions, resulting in a slight extension or compression 
of the Ti-0 bond pairs, and the "scissors" -type distor- 
tions of the Ti0 6 octahedron [HI QJ, QJ] . Th e Wyckoff 



positions of the atoms and their site symmetries in the or- 
thorhombic unit cell are listed in Table 1. The irreducible 
representations (in the Pbnm notation) corresponding to 
various atomic sites in the orthorhombic structure that 
follow from the character tables of the point groups are 
presented in the right column of Table 1. The total num- 
bers of the modes are grouped according to their op- 
tical activity. Among the total number of 60 T-point 
phonons, 24 (7A g +7Bi ff +5B 2£ ,+5B 3ff ) arc Raman-active 
modes, 25 (7Bi u +9B 2u +9B 3u ) are IR-activc modes, 8 
(A u ) are silent modes, and 3 (B lu +B 2u +B3 U ) are acous- 
tic modes. 



III. SAMPLES AND EXPERIMENTAL 
APPROACH 

A. Crystal growth and characterization 

It has to be emphasized that YTi03+,5 single crystals 
are always slightly off-stoichiometric, having an excess of 
oxygen above the exact stoichiometry formula. As to our 
knowledge, a minimal level record of S = 0.009 belongs to 
sing le crystals of YTi03+«5 grown by Czochralski method 
3l| . In this study, single crystals of YTi03 were grown 
by the floating zone method in a reducing atmosphere 
(Ar/H 2 = 50/50). More details are given in Ref. [H|. 
The oxygen excess in our samples was measured by the 
Differential Thermal Analysis and Thermal Gravimctry 
Analysis and estimated at a level less than S = 0.013. The 
same level of the oxygen off-stoichiometry 5 = 0.01 was 
determined in the single crystals of YTi03+<5, grown by 
the floating zone method in a fairly strong reducing con- 
dition (Ar/H 2 = 70/30) by Okimoto et al. [H. The op- 
tical conductivity spectra measured on our crystals and 
those grown by Okimoto et al. are very well agree with 
each other in a wide spectral range, from the phonon fre- 
quencies up to the deep UV frequencies (see Fig. 7 in Ref. 
[2l| and Fig. 2 in Ref. [H]). Therefore, we conclude that 
our single crystals of YTi03 grown by the floating zone 
method are of the same quality, verified by the optical 
probe. Moreover, the Ti 3+ ions in the YTi03 lattice are 
unstable in oxygen atmosphere, and fully oxidize to Ti 4+ 
above 300 °C in flowing oxygen (see Fig. 2 in Ref. [2ll]). 
This may lead to the aging effects in air, associated with 
the formation of the titanium oxides with the valence of 
Ti increased above 3+ and simultaneous incorporation 
of extra oxygen (interstitials). The aging effects and the 
oxygen off-stoichiometry in YTi03+a crystals have a deep 
impact on the optical conductivity spectra and their tem- 
perature dependencies [2l|, [H| . If the samples are kept in 
air for extending periods of time, the weak satellites can 
appear in the phonon modes, as arising from the coupling 
with the local oxygen modes associated with oxygen in- 
terstitials, as in the case of the ab surface, reported in 
Ref. [2l|. Therefore, an accurate study of the optical 
properties in YTi03+,5 crystals should be always referred 
to the oxygen off-stoichiometry 6. 
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The samples were further characterized by magnctom- 
etry using a superconducting quantum interference de- 
vice. For the almost stoichiometric YTi03+o.oi3 sample 
studied here, we estimate a Curie temperature of Tc = 
30 K, an inflection point in the temperature-dependent 
magnetization is observed at 27 K. The determined sat- 
urated magnetic moment in the easy direction at 7 T is 
fi c = 0.84 fiB, whereas fx b = 0.82 \xb in the hard direc- 
tion [2l[ . These our data are in good agreement with the 
magnetization data on the nearly oxygen stoichiometric 
YTi0 3 single crystals, which referred to T c of 29 K Q 
and 30 K [35|,[36[. 

At the same time x-ray diffraction measurements show 
that the YTiC-3 single crystals used for our experiments 
are untwinned, and posses an excellent crystal quality 
with mosaicity less than 0.03°. At room temperature, the 
lattice parameters are a = 5.331(3), b = 5.672(4), and c = 
7.602(6) A, slightly different from those of Czochralski- 
grown single crystals reported by Maclean et al. [3l1 |. 
The samples were aligned along the principal axes and 
cut in the form of a parallelepiped with dimensions ~ 
3x3x3 mm 3 . 



B. Ellipsometry approach 

For optical measurements the sample surfaces were 
polished to optical grade. The technique of ellipsometry 
provides significant advantages over conventional reflec- 
tion methods in that (i) it is self-normalizing and does 
not require reference measurements and (ii) ei(w) and 
£2(0;) are obtained directly without a Kramcrs-Kronig 
transformation [371 . l38l . [39j. For ellipsometry measure- 
ments at phonon frequencies 50-700 cm" 1 (0.006-0.087 
eV) we used the home-built ellipsometer in combina- 
tion with a fast-Fourier transform interferometer (Brukcr 
1FS 66v/S) at the infrared beam line of the ANKA syn- 
chrotron at the Forschungszentrum Karlsruhe, Germany. 
The high brilliance of the synchrotron light source due 
to the small beam divergence enabled us to perform very 
accurate ellipsomctric measurements in the far-IR range 
with incidence angles ranging from 65° to 85°. For low- 
temperature measurements the sample was mounted on 
the cold finger of a helium flow cryostat. With only 
a single angle of incidence, the raw experimental data 
are represented by real values of the ellipsometric angles, 
^ and A, for any wave number. These values are de- 
fined through the complex Fresncl reflection coefficients 
for light polarized parallel (r p ) and perpendicular (r s ) to 
the plane of incidence 

tan#e iA = ^. (1) 

To determine the complex dielectric response e a (ui), 
e&(w), and e c (uj) of the YTiC>3 crystal, we measured four 
high-symmetry orientations of the ac and be surfaces, 
with a or b axes aligned either parallel or perpendicu- 
lar to the plane of incidence of the light, respectively. 



Note here that we have used an overdetermined data set 
to check the consistency of our results. A nonlinear fit- 
ting procedure has been applied to extract point-by-point 
the complex dielectric response throughout the covered 
spectral range ■ In the following, we present the com- 
plex dielectric response e(u>) extracted from the raw el- 
lipsometry spectra, ^S>(uj) and A(w), according to this 
approach. 



IV. RESULTS AND DISCUSSION 

A. Infrared spectra of YTiC>3 

Figures 2-4 show the dielectric functions e(uj) and their 
inverses — l/e(w) for the three principal crystallographic 
axes of the YTiOs crystal, measured with an angle of 
incidence of 75° at T = 25 K. As discussed in Sec. II, a 
factor-group analysis for the orthorhombic crystal struc- 
ture of YTiOs (space group Pbnm, D^) with 4 f.u./unit 
cell yields a total number of 25 (7Bi u +9B2 M +9B3 U ) IR- 
active modes. Our polarized ellipsometry measurements 
allowed us to observe all symmetry-allowed IR-active 
phonon modes. One can easily identify six phonon modes 
in the polarization along the c axis (vibrations of Bi„ 
symmetry), eight phonon modes in the b axis (B2«), and 
eight phonon modes in the a axis (B3„). In addition, we 
were able to identify very weak features in our raw ellip- 
sometry spectra, \l/(w) and A(w), on the high-frequency 
side in each polarization (not shown), which we associate 
with Bi u (7), i?2«(9), and B3„(9) symmetry modes. Using 
the classical dispersion analysis we fit the complex dielec- 
tric functions e a {ui), £fc(w), and e c (u>) (and their inverses 
independently) with a set of Lorenzian oscillators 



T^—i ( 2 ) 

ujj — u) 1 — iwyj 



where u>j, jj, and Sj are the resonant frequency, width, 
and dimcnsionlcss oscillator strength of the jth oscilla- 
tor, and €00 is the core contribution from the dielectric 
function (4l| . The resulting parameters determine u>TOj 7 
7TOj, and Sj of the transverse TO phonon modes (and 
uJLOj and jhOj of the longitudinal LO phonon modes). 
We determine the anisotropic static dielectric constants 
eo from our fit at 50 cm -1 (these values are given in Table 

IV), and estimate the anisotropic high-frequency dielec- 

2 

trie constants defined by the relation e = e LlIm • 

In Table II we list all Lorentzian parameters, together 
with and e^, which show slightly different values, 
according to the accuracy of the present study. 
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B. Theoretical calculations 

1. Shell model approximation and potential parameters 

To determine the cigcnfrequcncics and eigenvectors of 
the normal modes and to make an assignment in the ex- 
perimental phonon spectra, we have performed lattice- 
dynamics calculations for the orthorhombic YT1O3 struc- 
ture in the Pbnm (D^h) space group in the framework of a 
shell model. The shell model approach has been success- 
fully applied earlier to study various physical properties 
and lattice dynamics of transition metal oxides, includ- 
ing high-T c superconductors [H, E3, S3 and LaMn0 3 
and YM11O3 [45I . liijl . ItH ]. Our calculations were per- 
formed in the GULP code [25|, using the lattice con- 
stants and fractional coordinates of the atoms reported 
for the Czochralski-grown single crystals by Maclean et 
al. [3l|. In the shell model approximation, the lattice 
is considered as an assembly of polarizable ions, repre- 
sented by massive point cores and massless shells cou- 
pled by isotropic harmonic forces defined by the spring 
constants. The interaction includes contributions from 
Coulomb and short-range interactions. The short-range 
potentials used for the shell-shell interactions are of the 
Buckingham form 

Vtj = Aijexpt-r/pij) - C^/r 6 . (3) 

The 0-0 (shell-shell) interactions are taken as those of 
typical oxides [H, EH • We adopt the formal ionic charges 
here. The forces acting on the atoms in the lattice were 
minimized by varying the Buckingham potential param- 
eters, spring constants, as well as the core and shell 
charges. The crystalline structure was allowed to relax 
to equilibrium conditions under the symmetry (Pbnm) 
restrictions, the cores and shells were allowed to move 
separately. The resulting Buckingham potential param- 
eters, shell charges (Y) and force constants (k) are pre- 
sented in Table III. The criterion for a successful fit was 
good agreement with the structural parameters [3l| , pro- 
vided that approximately zero forces are applied for all 
atoms, and with the experimental anisotropic static eo 
and high-frequency dielectric constants, determined 
in the present spectroscopic cllipsomctry study (as Table 
IV illustrates). 



2. Calculated Y -point IR- active modes in YTiOz 

A factor group analysis of the zone-center phonon 
modes in the orthorhombic YTi03 (Pbnm) structure 
yields 25 IR-active phonon modes, which have their pro- 
genitors in the Brillouin-zone of the cubic perovskite 
ABO3 (Pm3m) structure Among them, there are 

nine phonon modes, 3Bi„ + 3B2« + 3B3„, emanat- 
ing from the three IR-active Fi u cubic triplets. These 
are commonly referred to the external mode (in which 
the BOg octahedron vibrates against the A atoms), the 



bending mode (B-O-B bond-angle modulation), and the 
stretching mode (B-0 bond-length modulation). In ad- 
dition, two IR-active modes of B2« and B3„ symmetry 
and one silent mode of A„ symmetry stem from the cubic 
F 2u silent triplet. The rest orthorhombic IR modes arise 
from X, M, and R points of the cubic Brillouin-zone, 
and become active due to Brillouin-zone folding in the 
orthorhombic phase (50j . However, in this study we are 
not aimed at a detailed description of the orthorhombic 
phonon eigenvectors in terms of the corresponding cu- 
bic eigenvectors. Here we present results of our lattice- 
dynamics shell-model (SM) calculations and make an as- 
signment of the IR-active phonon modes wtoj observed 
in the polarized ellipsometry spectra e(oj) in the YTi03 
single crystal (see Figs. 2-4 and Table II). 

We tentatively divide the observed phonon modes into 
three groups centered at about the characteristic frequen- 
cies of external (170 cm -1 ), bending (340 cm -1 ), and 
stretching modes (540 cm -1 ) of the rare-earth perovskite 
oxides. Correspondingly, the eigenvector patterns for all 
normal modes of each symmetry Bi u , B2 U , and B3„, as 
derived from the SM lattice-dynamics calculations, are 
presented in Figs. 5-7. Here only major components of 
the eigenvectors (the whole list is given in Table V) along 
the principal crystallographic directions are displayed for 
20 atoms in the unit cell, of different symmetries and 
in all incquivalent positions. At each pattern the calcu- 
lated eigenfrcquencies are indicated. As one can notice 
from Table VI, the calculated cigcnfrequcncics generally 
well agree with the experimental TO phonon frequencies 
cuTOj- The character of the B iM , B 2u , and B 3u modes is 
determined according to the dominant contribution from 
the collective vibrations of Y, Ti, 01, or 02 atoms, polar- 
ized along the c-, b-, and a-axis, respectively. However, 
mixing due to out-of-phase vibrations of ions of equal 
symmetry in inequivalcnt positions often does not allow 
the determination of a unique character for a given set 
of equal symmetry modes. 

At low frequencies, near the characteristic frequency 
of the cubic perovskite external mode, one can identify 
three modes, B lu (l) : 167 cm -1 , B 2u (2) : 240 cm -1 , 
and B3 U (2) : 209 cm" 1 , which exhibit very high oscil- 
lator strengths in the experimental IR phonon spectra, 
in agreement with the SM calculations (see Table VI). 
The eigenvectors of the Bi„(l) mode show that the oscil- 
lator strength of this mode is determined by vibrations 
of the TiOe octahedra against the Y atoms along the c- 
axis. The eigenvectors of the B2„(2) and B3„(2) modes 
show that these normal modes involve dominant contri- 
bution from the collective Y-Ti vibrations, where Y and 
Ti atoms move in the opposite directions, polarized in 
accordance with the mode symmetries along the b- and 
a-axis, respectively. The eigenvectors of the B2 U (1) : 139 
cm -1 and B3„(l) : 154 cm -1 modes indicate that these 
modes are primarily due to out-of-phase vibrations of dif- 
ferent Y atoms and out-of-phase vibrations of different Ti 
atoms, polarized respectively along the a- and 6-axis, not 
coincident with the mode symmetries. In agreement with 
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the SM calculations, the observed phonon modes B2„(l) 
and B 3tl (l) exhibit very small oscillator strengths. An- 
other weak mode is the c-axis Bi„(2) : 204 cm' 1 , which 
is mainly due to 6-axis in-phase vibrations of Ti atoms in 
one plane and out-of-phase vibrations in different planes. 

At medium phonon frequencies near the cubic per- 
ovskitc bending mode, the strong c-axis polarized modes 
B lu (3) : 324 cm -1 and B lu (5) : 387 cm -1 are mainly 
due to Ti-02-Ti bending vibrations (the Bi„(3) mode 
has small admixture of 01 character due to the move- 
ment of 01 atoms in-phase with Ti atoms). The &-axis 
polarized modes B2„(3) : 273 cm -1 and B2 U (7) : 461 
cm -1 involve collective vibrations of 02 and Y atoms 
moving in the opposite directions. In turn, the strong 
a-axis polarized modes B 3u (3) : 316 cm' 1 and B 3tl (6) : 
426 cm -1 also involve collective vibrations of 02 atoms 
(the B 3tl (3) mode has in addition small contribution due 
to in-phase vibrations of 01 atoms and out-of-phase vi- 
brations of Y atoms). The character of the mode B2 U (4) 
: 308 cm" 1 is defined by the dominant contribution from 
01 vibrations, moving out-of-phase with 02 atoms. The 
modes B2„(5) : 364 cm -1 and B 3tl (5) : 344 cm -1 can 
be assigned as Ti-Ol-Ti bending modes. Other modes 
in this range B 2tl (6) : 380 cm -1 and B 3u (7) : 516 cm" 1 
are also determined by 01 collective vibrations, polarized 
according to the mode symmetries. We note here that 
the two 6-axis modes B2 U (5) : 364 cm -1 and B2 U (6) : 
380 cm -1 having very close frequencies could be reversed 
in the assignment, as it is explained below. The c-axis 
Bi„(4) : 349 cm -1 has very small oscillator strength as 
it can be assigned to the a-axis out-of-phase vibrations 
of Ti atoms. 

In the third group, the strong modes Bi„(6) : 545 
cm -1 , B 2u (8) : 528 cm" 1 , and B 3 „(8) : 554 cm -1 are ob- 
served in the experimental IR phonon spectra. Therefore, 
these IR-active phonon modes of orthorhombic YTi0 3 
should be most closely related to the cubic perovskite 
stretching modes. Indeed, the eigenvector patterns show 
that the Bi u (6) mode is due to the c-axis Ti-01 vibra- 
tions resulting in the Ti-01 bond-length modulation, and 
the oscillator strengths of the B2«(8) and B 3u (8) normal 
modes are determined by the b- and a-axis Ti-02 collec- 
tive out-of-phase vibrations. In addition, we were able to 
identify very weak features in our ellipsometry spectra 
on the high-frequency side in each polarization, which 
we associate with the Bi„(7) : 595 cm" 1 , £?2«(9) : 577 
cm" 1 , and B 3u (9) : 576 cm -1 symmetry modes. We 
note that for the Bi„(6) and Bi„(7) modes there is a 
good agreement between the frequencies and oscillator 
strengths, predicted from the calculations and those ob- 
served (see Table VI). The Bi„(7) normal mode is weak, 
because it is mostly due to the different symmetric ab- 
polarized out-of-phase oscillations of Ti and 02 atoms. 
However, our SM calculations predict that the strong in- 
plane modes of B2 U and B% u symmetry appear at the 
high-frequency boundary of the phonon spectra. In this 
case, in our assignment we relate the modes of similar 
oscillator strengths. In agreement with our SM calcula- 



tions the B2„(9) and B 3tl (9) normal modes are largely 
due to the symmetric a- and 6-polarized out-of-phase vi- 
brations of Ti and 01 atoms, respectively. Regarding the 
symmetry of these modes, their oscillator strengths is ex- 
pected to be small, in consistency with our experimental 
observations. 



3. Comparison with first-principles results 

In order to independently support the assignment of 
the IR-active phonon modes in YTi0 3 , we have calcu- 
lated the T-point normal modes by means of the ab initio 
approach within the density functional theory using the 
VASP software package [27], [H| . A detailed description 
of these calculations will be given elsewhere [26| . Here we 
only mention the eigenfrequencies of the normal modes 
resulting from the ab initio calculations: a-axis 9 B 3tl 
normal modes are calculated at the frequencies 137, 175, 
298, 312, 330, 379, 391, 473, and 519 cm" 1 ; 6-axis 9 B 2u 
normal modes at 125, 211, 250, 291, 337, 359, 440, 485, 
and 500 cm" 1 ; and c-axis 7 Bi u normal modes at 144, 
182, 289, 336, 352, 463, and 517 cm" 1 . In general, the fre- 
quencies and eigenvectors of the normal modes predicted 
by the ab initio lattice-dynamics calculations are in a 
reasonable agreement with the results of the shell model 
lattice-dynamics calculations reported in detail here, es- 
pecially for the low-frequency and high-frequency modes, 
near the characteristic frequencies of the cubic perovskite 
external and stretching modes, respectively. The main 
characters of the modes near the characteristic cubic per- 
ovskite bending mode (where the phonon frequencies are 
densely spaced and therefore the results are very sensi- 
tive to details) are also in satisfactory agreement with the 
shell model lattice-dynamics calculations. For example, 
the assignments of the 6-axis modes B2«(3) : 273 cm" 1 
and B2„(4) : 308 cm -1 , the a-axis mode B 3u (6) : 426 
cm" 1 and the c-axis mode Bi u (3) : 324 cm -1 are well 
reproduced. 

4- Effective charges in YTiOz 

The static effective charges Z^e related to the bonding 
ionicity can be defined for a ternary compound as follows 

Y\ut Qj - w%o S ] = y ^(Z^/M,, (4) 

where j denotes the phonon mode with TO frequency 
i^TOj and LO frequency wloj , V c is the unit cell volume, 
Mi corresponds to the mass of the ith atom. The static 
effective charges of oxygen Zoe can be estimated rather 
accurately from this relation, as the contribution of the 
terms related to oxygen is dominant on the right hand 
side of Eq. (4) due to the small mass of oxygen as com- 
pared with the masses of other constituents. Using the 
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frequencies wtoj and uJhOj determined from the ellip- 
sometry measurements listed in Table II, we have found 
that the static effective charges of oxygen in orthorhom- 
bic YTi03 are nearly isotropic: Zo« a = -1.14 | e |, Zo« b 
= -1.17 | e |, and Zo« c = -1.20 | e |. These values are 
typical for the entire family of perovskite oxides. 

The Born effective charge of an ion represents the dy- 
namical contribution to the effective charge. It is con- 
trolled by long-range Coulomb interactions. The dynam- 
ical effective charge is represented by a tensor, since it is 
defined by the change of macroscopic polarization created 
in direction /3, induced by the periodic displacement r K a 
of the sublattice of atoms k in direction a at the linear 
order, times the unit cell volume 

Q K , a /3 = Vcjr — • (5) 

For the case of ABO3 compounds, the values of the Born 
effective charges deviate substantially from the static ef- 
fective charges. The ah initio calculations report anoma- 
lously large values of dynamical effective charges for 
Ti and O ions in ferroelectric perovskites KNbC>3 and 
BaTi0 3 HI, [H. The charges of Ba and Ti atoms in 
the orthorhombic structure Q% a and of BaTiC>3 are 
nearly isotropic and vary, respectively, in the range from 
+2.72 to +2.77 and from +5.59 to +6.8. In turn, for 
oxygen, the effective charges exhibit predominantly uni- 
axially anisotropic character, where the values of Qq 
vary from -1.91 to -2.04 and Qq ± vary from -4.89 to - 
5.45, and refer, respectively, to a displacement of oxygen 
in the a6-plane and along the c-axis [54] ] . 

The oscillator strengths Sj of IR-active phonon modes 
are related to the dipole moment arising from the dis- 
placements of ions involved in a vibration 

j V c An, V. M iU l ■ M,i 

Here Ql corresponds to the dynamical effective charge 
of the ith atom vibrating in j mode with the TO fre- 
quency wxoj: an d Uij represents the displacements re- 
lated to the eigenvector Xj by ity = Xij/VMi. In Table 
VI we present calculated squares of net dipole moments 
E< Qi u ij] 2 arising from the displacements of all ions in- 
volved in the different symmetry IR-active vibrations, us- 
ing the formal ionic charges that were adopted in our 
shell model. These values are juxtaposed with the corre- 
sponding values estimated from the experimental phonon 
parameters of u>TOj an d oscillator strengths Sj given in 
Table II. As one can see, the qualitative agreement be- 
tween the model and the experiment is quite good in 
the sense that the model correctly predicts the strongest 
modes in the phonon spectra and the weak modes are 
also predicted to be weak. However, some of the modes, 
for example, two strong 6-axis modes B2 U (3) : 273 cm -1 
and B 2 „(6) : 380 cm -1 , two strong a-axis modes B 3u (2) 
: 209 cm -1 , B 3 „(6) : 426 cm -1 , and one weak mode 
63^(7) : 516 cm -1 , and the strong c-axis mode Bi u (3) 



: 324 cm -1 show largest discrepancy between the exper- 
imental values and calculated in the framework of the 
adopted shell model. 

An improved description of the experimental oscilla- 
tor strengths of all phonon modes compared to those 
with formal ionic charges can be achieved by introducing 
Born effective charges. These effective charges can be ex- 
tracted by inversion the linear equations Eq. (6) by mak- 
ing use of the eigenvectors, providing the experimental 
phonon parameters of WxOj & n d Sj are available. We note 
here that for each polarization together with the electric 
charge neutrality condition, we have an over determined 
data set for the four unknown effective charges. However, 
the inversion of Eq. (6) results in large errors if it is ap- 
plied for the modes which have eigenvector components 
close to zero, in particular, for the weak modes (see Ta- 
bles V and VI), as the determinant of the matrix is then 
close to zero. We apply Eq. (6) to the three strong modes 
in each polarization and determine the corresponding dy- 
namical effective charges, which are listed in Table VII. 
In reasonable agreement with the effective charges of Ti 
atoms in the orthorhombic structure of BaTiC>3, with 
the nominal valence of Ti atoms +4 [54[ , we have found 
that in YTiC>3 the dynamical effective charges of are 
nearly isotropic, varying in the range from +3.9 to +5.00. 
At the same time, the effective charges of Y and 02 atoms 
exhibit a-axis anisotropy: Qyii a =+1.5 and ,i=-1.9, 
compared to the corresponding effective charges Qy|| bc 
of +3.4^- +3.9andQ^ i2!|fc of-2.5^-2.7. The dynamical 
effective charges of oxygen Qq i ± , related to a displace- 
ment along the c-axis, vary in the range from -2.3 to -2.8 
in YTi03 and do not show anomalously large values like 
in ferroelectric BaTi03 [54]. In Table VI the calculated 
values of the Q*Uij} 2 are juxtaposed with the corre- 
sponding experimental values and those obtained from 
the calculations with formal ionic charges. One can no- 
tice that, in general, the agreement with the experimental 
values is improved (except for some of the weak modes, 
which are less accurately predicted). The description 
with the effective charges give also an indication that two 
6-axis modes B 2m (5) : 364 cm -1 and B 2m (6) : 380 cm -1 
having very close frequencies should be reversed in the 
assignment. The agreement for the a-axis modes B3„(6) 
: 426 cm" 1 and B3„(7) : 516 cm" 1 is not substantially 
improved by introducing Born effective charges. A pos- 
sible reason could be a strong dependence of the mode 
B3«(7) on strain in the real crystal structure. As a result, 
most of its oscillator strength might be transferred to the 
strongest a-axis mode B3„(6). We note that the 6-axis 
mode B2 M (3) : 273 cm" 1 is less successfully predicted, 
however, the character of this mode is well reproduced 
by the present ab initio calculations ;26;]. 

V. CONCLUSIONS 

In summary, we present optical phonon spectra of the 
orthorhombic YTi03 single crystal as obtained by spec- 
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troscopic cllipsometry and suggest the assignments of all 
symmetry-allowed zone-center IR-active phonon modes 
by comparing them with the results of the shell model 
calculations. We note that the shell model elaborated 
in this work may possess certain advantages, as the pa- 
rameters of the model were fitted not only to the or- 
thorhombic crystal structure, but also to the experimen- 
tal anisotropic static £q and high-frequency eoo dielectric 
constants, determined in the present cllipsometry study. 
Consequently, the calculated dipole moments (or oscilla- 
tor strengths) associated with displacements of the atoms 
involved in the IR-active normal modes in YTi03 are in 
good agreement with the experimental data. The normal 
modes predicted in the framework of the present shell 
model are in good agreement with the ab initio lattice- 
dynamics calculations [261 ] - Thus, our study gives us a 
coherent picture of the lattice vibrations in YTiC>3, con- 
sistent with both experimental data and with full elec- 
tronic structure calculations. We have not considered the 
influence of magnetic interactions on the lattice dynam- 
ics, but as the ferromagnetic coupling in YTi03 is weak 
it should not affect our results substantially. 

In order to describe the lattice dynamics of complex 
oxides like the rare-earth titanates, one can easily imag- 
ine a need to go beyond the shell model combined with 
interatomic potentials. However, the shell model itself 
is a general model within the dipole and harmonic ap- 
proximations [55l | , and there have been wide-ran ging and 
successful calculations for simple binary oxides [561 ] . in- 
cluding predictions of surface structure and vibrations, 
off-center ions, charge transfer transitions, and defect 
equilibria for non-stoichiomctric oxides. What was un- 
expected is the extent to which the oxygen-oxygen po- 
tentials seem to be transferrable to the more complex 
oxides. One might have expected the Jahn- Teller effect 
or charge transfers and lack of ionicity to lead to radi- 
cal changes. This seems not to be the case: we find the 
0-0 potentials based on MgO and similar binary oxides 
essentially the same for LaMn03, for Bi-based oxide su- 
perconductors, and for YTi03. Some of the reasons for 
this are discussed in relation to oxide superconductors in 
Ref. [57j] ■ As a result, the shell model combined with 
interatomic potentials continues to offer a powerful route 
to understanding experimental data for complex oxides. 



The comprehensive study described here provides a 
solid foundation for further quantitative exploration of 
the interplay between spin, orbital and lattice degrees of 
freedom in the titanates. We found that in YTi03 the 
effective charges of Ti are nearly isotropic and compa- 
rable to the effective charges of Ti in other perovskite 
compounds. At the same time the effective charges of 
Y and 02 atoms exhibit a-axis anisotropy: Qyii a = +l-5 
and Qo 2 ii a = -l-9, compared to the corresponding effective 
charges 2 Q^ ||bc (+3.4 - +3.9) and Q*^^ (-2.5 - -2.7). 
As the effective charges are a sensitive tool for analyz- 
ing dynamic changes of orbital hybridization, we suggest 
strong hybridization effects between Y and et6-plane oxy- 
gen along the a-axis. Hybridization of Y 5d and O 2p 
states has indeed been predicted by density functional 
calculations [l3|. Then, according to our assignment, the 
modes B3„(3) : 316 cm -1 and B3„(6) : 426 cm -1 , involv- 
ing 02 and Y displacements along the a axis, will proba- 
bly exhibit anomalous temperature dependencies. Some 
of the Raman- active modes related to the "scissors" -type 
distortions of the TiOf, octahedra, which are especially 
important in YTi03 [14| . might also exhibit anomalous 
temperature dependence in the vicinity of Tc- As in 
the real structure of YTi03 these distortions are most 
pronounced in the 6c-plane, we expect that also some of 
the IR-active modes involving 01 oxygen displacements 
along the b axis and 02 oxygen displacements along the 
c axis will exhibit anomalous behavior around the fer- 
romagnetic transition temperature. According to our as- 
signment, the possible candidates are B2 M (4) : 308 cm -1 , 
B 2u (6) : 380 cm" 1 , Bi u (3) : 324 cm" 1 , and B lu (5): 387 
cm" 1 . 



VI. ACKNOWLEDGMENTS 

The authors thank J. Gale for making available Gen- 
eral Utility Lattice Program (GULP) used in the present 
calculations. We are grateful to R. Evarestov, E. Ko- 
tomin and O. Dolgov for fruitful discussions. We ac- 
knowledge Y.-L. Mathis at ANKA for his support during 
the cllipsometry measurements. 



[1] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998). 
[2] T. Mizokawa, D. I. Khomskii, and G. A. Sawatzky, Phys. 

Rev. B 60, 7309 (1999). 
[3] H. Sawada and K. Terakura, Phys. Rev. B 58, 6831 

(1998). 

[4] B. Keimer, D. Casa, A. Ivanov, J. W. Lynn, M. v. 

Zimmermann, J. P. Hill, D. Gibbs, Y. Taguchi, and Y. 

Tokura, Phys. Rev. Lett. 85, 3946 (2000). 
[5] G. Khaliullin and S. Maekawa, Phys. Rev. Lett. 85, 3950 

(2000). 



[6] G. Khaliullin and S. Okamoto, Phys. Rev. Lett. 89, 
167201 (2002). 

[7] C. Ulrich, G. Khaliullin, S. Okamoto, M. Reehuis, A. 

Ivanov, H. He, Y. Taguchi, Y. Tokura, and B. Keimer, 

Phys. Rev. lett. 89, 167202 (2002). 
[8] J. Akimitsu, H. Ichikawa, N. Eguchi, T. Miyano, M. 

Nishi, and K. Kakurai, J. Phys. Soc. Jpn. 70, 3475 

(2001). 

[9] M. Mochizuki and M. Imada, J. Phys. Soc. Jpn. 70, 2872 
(2001); M. Mochizuki, J. Phys. Soc. Jpn. 71, 2039 (2002); 
M. Mochizuki and M. Imada, New J. Phys. 6, 154 (2004). 



8 



H. D. Zhou and J. B. Goodenough, J. Phys.: Condens. 
Matter 17, 7395 (2005). 

M. Cwik, T. Lorenz, J. Baier, R. Miiller, G. Andre, F. 
Bouree, F. Lichtenberg, A. Freimuth, R. Schmitz, E. 
Muller-Hartmann, and M. Braden, Phys. Rev. B 68, 
060401(R) (2003). 

I. V. Solovyev, Phys. Rev. B 69, 134403 (2004); I. V. 
Solovyev, Phys. Rev. B 74, 054412 (2006). 

E. Pavarini, S. Biermann, A. Poteryaev, A. I. Lichten- 
stein, A. Georges, and O. K. Andersen, Phys. Rev. Lett. 
92, 176403 (2004); E.Pavarini, A. Yamasaki, J. Nuss, and 
O. K. Andersen, New J. Phys. 7, 188 (2005). 
A. A. Mozhegorov, A. V. Larin, and A. E. Nikiforov, 
Phys. Met. Metallogr. 105, 219 (2008). 
R. Schmitz, O. Entin-Wohlman, A. Aharony, A. B. Har- 
ris, and E. Miiller-Hartmann, Phys. Rev. B 71, 214438 
(2005). 

J. E. Greedan, J. Less-Common Met. Ill, 335 (1985). 
A. C. Komarek, H.Roth, M Cwik, W.-D. Stein, J. Baier, 
M. Kriener, F. Bouree, T. Lorenz, and M. Braden, Phys. 
Rev. B 75, 224402 (2007). 

C. Ulrich, A. Gossling, M. Griininger, M. Guennou, 
H. Roth, M. Cwik, T. Lorenz, G. Khaliullin, and B. 
Keimer, Phys. Rev. Lett. 97, 157401 (2006). 

V. Fritsch, J. Hemberger, M.V. Eremin, H.A. Krug von 
Nidda, F. Lichtenberg, R. Wehn, and A Loidl, Phys. Rev. 
B 65, 212405 (2002). 

J.-G. Cheng, Y. Sui, J.-S. Zhou, J.B. Goodenough, and 
W.H. Su, Phys. Rev. Lett. 101, 087205 (2008). 
N. N. Kovaleva, A. V. Boris, P. Yordanov, A. Maljuk, E. 
Bruchner, J. Strempfer, M. Konuma, I. Zegkinoglou, C. 
Bernhard, A. M. Stoneham, and B. Keimer, Phys. Rev. 
B 76, 155125 (2007). 

W. Knafo, C. Meingast, H. v. Lohneysen, A. V. Boris, P. 
Popovich, N. N. Kovaleva, P. Yordanov, A. Maljuk, R. 
K. Kremer, and B. Keimer, submitted to Phys. Rev. B. 
J. Zhang, P. Dai, J. A. Fernandez-Baca, E. W. Plummer, 
Y. Tomioka, and Y. Tokura, Phys. Rev. Lett. 86, 3823 
(2001). 

E. Saitoh, Y. Tomioka, T. Kimura, and Y. Tokura, J. 
Phys. Soc. Jpn. 69, 2403 (2000). 

J. D. Gale, Phil. Mag. B 73, 3 (1996); J. Chem. Soc. 
Faraday Trans. 93, 629 (1997). 
L. Capogna, unpublished. 

G. Kresse and J. Furthmuller, Phys. Rev B 54, 11169 
(1996). 

K. Parlinski, in Neutrons and Numerical Methods, AIP 
Conference proceedins, 479 (1999). 

International Tables for X-ray Crystallography, Ed. by 
Norman F. M. Henry and Kathleen Lonsdale (Kynoch 
Press, Birmingham, 1969). 

A. M. Glazer, Acta Crystallogr., Sect. B: Struct. Crys- 
tallogr. Cryst. Chem. B28, 3384 (1972). 

D. A. MacLean, Hok-Nam Ng, and J. E. Greedan, J. 
Solid State Chem. 30, 35 (1979). 

Y. Okimoto, T. Katsufuji, Y. Okada, T. Arima, and Y. 

Tokura, Phys. Rev. B 51, 9581 (1995). 

A. G0"ssling, R. Schmitz, H. Roth, M. W. Haverkort, 

T. Lorenz, J. A. Mydosh, E. M"ller-Hartman, and M. 

Gr"ninger, Phys. Rev. B 78, 075122 (2008). 

J. D. Garret and J. E Greedan: Inorg. Chem. 20, 1025 

(1981). 

M. Tsubota, F. Iga, T. Takabatake, N. Kikugawa, T. 



[36] 

[37] 
[38] 
[39] 
[40] 

[41] 



Suzuki, I. Oguro, H. Kawanaka, and H. Bando, Physica 
B 281-282, 622 (2000). 

N. Tsuji, M. Ito, H. Sakurai, K. Suzuki, K. Tanaka, K. 
Kitani, H. Adachi, H. Kawata, A. Koizumi, H. Nakao, 
Y. Murakami, Y. Taguchi, and Y. Tokura, J. Phys. Soc. 
Japan 77, 023705 (2008). 

R. Henn, C. Bernhard, A. Wittlin, M. Cardona, S. 
Uchida, Thin Solid Films 313-314, 643 (1998). 
C. Bernhard, J. Humlicek, and B. Keimer, Thin Solid 
Films 455-456, 143 (2004). 

Handbook of Ellipsometry, Ed. by H. G. Tompkins and 

E. A. Irene (Springer, Berlin, 2005). 

J. Humhcek, R. Henn, and M. Cardona, Phys. Rev. B 

61, 14554 (2000); J. Humhchek, Philos. Mag. B 70, 699 

(1994). 

To extract information about the phonons we fitted the 
ellipsometry spectra in Figs. 2-4 using a simple Lorentz 
model for noninteracting harmonic oscillators. However, 
as one can see from Figs. 2-4, several phonon bands show 
weak asymmetry, with flatter rise and steeper fall of the 
imaginary part of e(u>). At the same time the —l/i(u>) 
spectra were found to be fitted fairly well with a Lorentz 
model. The asymmetry of optical phonon modes suggests 
mode coupling and interaction damping, which can be ac- 
counted for by a phenomenological model for coupled os- 

2 _ - 

cillators, e(w) = + E, Sj „, ilJ^L.- > with the asym- 



[42] 

[43] 
[44] 

[45] 



[46] 
[47] 



[49] 
[50] 

[51] 
[52] 



[53] 

[54] 

[55] 
[56] 

[57] 



— 'J J u>.— uj z — zuj~fj 

metry parameters <7j [4(J. We estimate that the level of 
asymmetry of all TO phonons in YTi03 is fairly weak 
and within the experimental uncertainties in our case. 
J. Prade, A. D. Kulkarni, F. W. de Wette, U. Schroder, 
W. Kress, Phys. Rev. B 39, 2771 (1989); A. D. Kulkarni, 
F. W. de Wette, J. Prade, U. Schroder, W. Kress, Phys. 
Rev. B 41, 6409 (1990). 

M. S. Islam, M. Leslie, S. M. Tomlinson, and C. R. A. 

Catlow, J. Phys. C 21, L109 (1988). 

N. N. Kovaleva, A. V. Boris, T. Holden, C. Ulrich, B. 

Liang, C. T. Lin, B. Keimer, C. Bernhard, J. L. Tallon, D. 

Munzar, and A. M. Stoneham, Phys. Rev. B 69, 054511 

(2004). 

M. N. Iliev, M. V. Abrashev, H.-G. Lee, V. N. Popov, Y. 
Y. Sun, C. Thomsen, R. L. Meng, and C. W. Chu, Phys. 
Rev. B 57, 2872 (1998). 
I. S. Smirnova, Physica B 262, 247 (1999). 

N. N. Kovaleva, J. L. Gavartin, A. L. Shluger, A. V. 
Boris, and A. M. Stoneham, JETP 94, 178 (2002). 
C. R. A. Catlow, W. C. Mackrodt, M. J. Norgett, A. M. 
Stoneham, Philos. Mag. 35, 177 (1977). 
V. N. Popov, J. Phys.: Condens. Matter 7, 625 (1995). 
Michel Couzi and Pham Van Huong, J. Chim. Phys. 69, 
1339 (1972). 

J. F. Scott, Phys. Rev. B 4, 1360 (1971). 
S. Tajima, T. Ido, S. Ishibashi, T. Itoh, H. Eisaki, Y. 
Mizuo, T. Arima, H. Takagi, and S. Uchida, Phys. Rev. 
B 43, 10496 (1991). 

R. Resta, M. Posternak, and A. Baldereschi, Phys. Rev. 
Lett. 70, 1010 (1993). 

Ph. Ghosez, X. Gonze, Ph. Lambin, J.-P. Michenaud, 

Phys. Rev. B 51, 6765 (1995). 

W. Cochran, Crit. Rev. Sol. St. Sci. 2, 1 (1971). 

A. M. Stoneham and J. H. Harding, Ann. Rev. Phys. 

Chem. 37, 53 (1986). 

A. M. Stoneham and L. W. Smith, J. Phys. Cond. Mat. 
3, 225 (1991). 




FIG. 1: Crystal structure of cubic perovskite oxides ABO3 
(space group Pm3m) where A is a rare-earth atom (yel- 
low sphere) and B is a transition-metal atom (red spheres) 
(left), and related orthorhombic structure of YTiOs (space 
group Pbnm), with ferromagnetic alignment of Ti spins shown 
schematically (right). 
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TABLE I: Wyckoff notations, atomic site symmetries, and irreducible representations for the atoms in orthorhombic Pbnm 



(Dft) YTi0 3 . 


Atom 


Wyckoff 
notation 


Site 

symmetry 


Irreducible representation 


Y 
Ti 
01 
02 


4(c) 
4(b) 
4(c) 
8(d) 


c" 

Qxy 

Ci 


2 A g +A u + 2Bi 9 + Bi„ + B2 S + 2B2 U + B3 9 + 2B3„ 

3 A u + 3Bi u + 3B2u + 3B3 U 

2A g + A u + 2Bi g + Bi u + B2g + 2B2 U + B3 S + 2B3„ 
3A g + 3A U + 3Bi 9 + 3Bi„ + 3B2 S + 3B2 U + 3B 3s + 3B3„ 



7A S + 8A„ + 7Bi 9 + 8Bi„ + 5B 2s + 10B 2u + 5B 3s + 10B 3u 
IYr = 7Bi u + 9B2 U + 9B 3u 

^Raman = 7A g + 7Bl 9 + 5B2 9 + 5B3g 
^silent = 8A U 

J- acoustic = Bl„ + B211 + B3 U 



TABLE II: Results of the dispersion analysis based on the three-parameter Lorenz model: frequencies and damping parameters 
lotO] and 7tOj (in cm" 1 ), and dimensionless oscillator strengths Sj of the TO phonon modes in the the polarized e(u>) spectra 
of the YTiOs single crystal (T = 25 K), eoo is the core contribution; frequencies and damping parameters cjlOj and "/LOj (in 
cm -1 ) of the LO phonon modes in the the -l/e(co) spectra. Definition of is explained in the text. 



b axis 





7toj 


Sj 




7LOj 


OJTOj 


7TOj 


Sj 


^LOj 


7LOj 




7tOj 


Si 




7LOj 


154 


5.5 


0.38 


155 


4.3 


139 


5.5 


0.53 


141 


4.0 


167 


3.8 


7.15 


196 


2.9 


209 


5.1 


3.77 


236 


5.2 


240 


8.8 


2.36 


253 


5.0 


204 


4.9 


0.30 


216 


6.5 


316 


6.9 


2.05 


331 


3.6 


273 


7.3 


0.60 


279 


9.8 


324 


8.6 


3.02 


346 


4.8 


335 


3.1 


0.12 


342 


8.5 


308 


5.3 


0.51 


313 


4.6 


349 


6.2 


0.08 


354 


9.5 


344 


9.7 


0.09 


351 


6.6 


364° 


10.0 


1.37 


369 


7.9 


387 


9.0 


2.37 


484 


7.6 


426 


7.7 


2.78 


508 


13.2 


380 


9.9 


2.59 


438 


6.9 


545 


7.2 


0.85 


695 


7.3 


516 


10.4 


0.07 


525 


7.5 


461 


8.8 


0.61 


496 


9.2 


595 


12.0 


0.01 


595 


12.5 


554 


7.1 


0.61 


683 


7.1 


528 


12.5 


0.82 


689 


10.6 












576 


11.3 


0.02 


576 


11.2 


578 


11.2 


0.02 


577 


9.6 













Coo — 3.5 


£00=4.0 


£oo=3.1 


£^=4.0 


e* —4 1 

c oo 


£00=3.8 



"the splitting of this TO phonon in b-axis is probably due to a 
cross talk of the a-axis with strong LO mode at 351 cm -1 . 



TABLE III: Buckingham potential parameters for the shell-shell and core-shell interactions, the ionic shell charges (Y) (in units 
of a free electron charge) and force constants (k) in orthorhombic YTiOs, r cu toff = 19 A. 



A(eV) p(A) C(eV-A- 6 ) 



O-O 




22764.0 




0.1490 




20.37 


Ti-Ol 




2902.8 




0.2837 




00.00 


Ti-02 




1240.1 




0.3313 




00.00 


Y-Ol 




3860.6 




0.2944 




00.00 


Y-02 




12769.1 




0.2495 




00.00 




Ion 




Y(|e|) 




k (eV -A 






oi 2 ~ 




-2.000342 




34.4 






02 2 " 




-3.271589 




35.7 
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TABLE IV: Experimental and calculated properties of YTi03. The lattice constants and fractional atomic coordinates present 
the data for the Czochralski- grown single crystals at T = 293 K [3lj . 

Properties Experiment Calculation 

a^A 5.316(2) 5.3124 

b, A 5.679(2) 5.6880 

c, A 7.611(3) 7.6072 
V, A 3 229.799 229.864 

x/a, y/b, z/c 

Y 0.9793(1),0.0729(1),0.25 0.9778,0.0692,0.25 

Ti 0.5,0.0,0.0 0.5,0.0,0.0 

01 0.121(1),0.458(1),0.25 0.1185,0.4485,0.25 

02 0.691(1),0.310(1),0.058(1) 0.6967,0.3055,0.0595 

Dielectric constants 

e£ 13?7 13.34 

e§ 2 13.6 14.02 

eg 3 17.7 17.52 

3.5-4.0 3.84 

eg 4.0-4.1 4.31 

e§| 3.1-3.8 3.21 



TABLE V: Calculated eigenfrequencies Uxc-j anc ^ components of eigenvectors Xij of the IR normal modes in orthorhombic 
YTi03. Signs of the components are referred to the atoms with fractional coordinates in the unit cell: Y(0. 9793, 0.0729, 0.25), 
Ti(0.5,0,0), 01(0.121,0.458,0.25), and 02(0.691,0.310,0.058). The eigenvector components used in the inversion of Eq. (6) are 
underlined. 
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TABLE VI: Suggested assignments of the TO phonon modes in the polarized ellipsometry spectra, WTOj, with the normal 
modes resulting from the shell-model lattice-dynamics calculations, ui^f. Squares of net dipole moments calculated with 
the eigenvector components the formal ionic charges Qi and the effective dynamical charges Q* are compared with the 
corresponding values estimated from the experimental phonon parameters tJTOj an d Sj . 
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TABLE VII: The dynamical effective charges of the atoms (in units of a free electron charge) calculated from experimental 
phonon parameters, LJTOj and Sj, using the eigenvector components Xij listed in Table V. Uncertainties of the effective charges 
are estimated by varying the experimental oscillator strength by ±10%. 
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FIG. 2: Real (red open circles) and imaginary (black open 
circles) parts for the measured a-axis complex dielectric re- 
sponse e(w) and its inverse -l/e(u>). The solid curves present 
the result of the fitting with a set of Lorenzian oscillators, 
which resonant frequencies lutoj and ajLOj are indicated (see 
also Table II). 




FIG. 3: The same as in Fig. 2 for the &-axis complex dielectric 
response e(w) and its inverse -l/e(w). 
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FIG. 4: The same as in Fig. 2 for the c-axis complex dielectric 
response e(u>) and its inverse -l/e(w). 




FIG. 5: Eigenfrequencies and eigenvector components of 
the normal modes of Bi u , B2u, and Bz u symmetries in or- 
thorhombic YTiOa at low frequencies, near the characteristic 
frequency of the perovskite external mode. 
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FIG. 6: Eigenfrequencies and eigenvector components of 
the normal modes of Bi u , B2u, and Bz u symmetries in or- 
thorhombic YTIO3 near the characteristic frequency of the 
perovskite bending mode. 




FIG. 7: Eigenfrequencies and eigenvector components of 
the normal modes of B\ u , B2 U , and B^ u symmetries in or- 
thorhombic YTi03 at high frequencies, near the characteristic 
frequency of the perovskite stretching mode. 



